hep-th/0207107 
CTP-MIT-3283 



Dynamics with Infinitely Many Time Derivatives 
and Rolling Tachyons 

Nicolas Moeller and Barton Zwiebach 



Center for Theoretical Physics 
Massachusetts Institute of Technology, 
Cambridge, MA 02139, USA 

E-mail: moeller, zwiebacli@mitlns.mit.edu 



Abstract 

Both in string field theory and in p-adic string theory the equations of motion involve 
infinite number of time derivatives. We argue that the initial value problem is qualitatively 
different from that obtained in the limit of many time derivatives in that the space of initial 
conditions becomes strongly constrained. We calculate the energy-momentum tensor and 
study in detail time dependent solutions representing tachyons rolling on the p-adic string 
theory potentials. For even potentials we find surprising small oscillations at the tachyon 
vacuum. These are not conventional physical states but rather anharmonic oscillations 
with a nontrivial frequency-amplitude relation. When the potentials are not even, small 
oscillatory solutions around the bottom must grow in amplitude without a bound. Open 
string field theory resembles this latter case, the tachyon rolls to the bottom and ever 
growing oscillations ensue. We discuss the significance of these results for the issues of 
emerging closed strings and tachyon matter. 
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1 Introduction and Summary 

Some of the most intriguing and fascinating aspects of string theory center around the 
role of locality. More precisely, string theory appears to be a consistent theory that is not 
local in the sense of local quantum field theory. Indeed, the covariant string field theory 
action governing the dynamics of the infinite set of spacetime fields contains spacetime 
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derivatives of all orders. Not only we have spatial derivatives of all orders, as in non- 
commutative field theory, but we also have time derivatives of all orders, and mixed 
derivatives of all orders. 

The possible difficulties with theories having a high number of time derivatives are 
familiar. They include possible violations of unitarity and causality, the appearance of 
spurious solutions, complications setting up an initial value problem, and difficulties quan- 
tizing the theory and/or finding a stable Hamiltonian [|l|. There is, however, evidence 
that in the case of string theory the higher derivative structure does not threat the consis- 
tency of the theory. The perturbative S-matrix of string theory is unitary. Even though 
covariant string field theory cannot be quantized using a Hamiltonian, path integral quan- 
tization is possible and elegantly done with the use of the Batalin-Vilkovisky formalism 
[|], ^, ^, 0. While infinitely many spatial derivatives are possible in non-commutative field 
theory, infinitely many mixed spatial and time derivatives appear to require string theory 
for consistency p, 0. Finally, the higher derivative structure of string theory is certainly 
very special - when passing to the light-cone gauge the theory becomes local in light-cone 
time. Indeed, light-cone string field theory has only first order time derivatives M. 



Motivated by the recent work of Sen on tachyon matter |T^, |Tl| we examine the 
rolling of the tachyon down its potential. Tachyon matter is a classical open string theory 
solution where at large times the tachyon T{t) approaches the tachyon vacuum T = oo 
with constant velocity T{t) . This solution represents a pressureless gas and its possible 
cosmological relevance has been studied in |T^. The evidence for such solutions comes 
mostly from the conformal field theory approach. Recently, related solutions have been 
obtained from considerations of boundary string field theory [l^, |14|. These actions. 



truncated to include all powers of first derivatives of fields exhibit similar, though not 
identical behavior. The purpose of the present paper is to examine tachyon dynamics in 
the setup of open string field theory [jl5l , where the equations of motion contain infinitely 
many time derivatives. Most of our work, however, will concentrate in the case of p- 
adic string theory This string theory is essentially a theory of a scalar field with 



infinitely many spacetime derivatives. Indeed many of the properties of the tachyon field 
of open string theory are also exhibited by the p-adic string model [|17|, |18[. One has 
a tachyon, and a potential with a locally stable minimum where the scalar field has no 
conventional excitations. Lump-like solutions exist representing p-adic D-branes of various 
dimensionalities. We will find here interesting phenomena regarding tachyon solutions in 
p-adic string theory, and we will also see that some solutions in string field theory are 
quite analogous to p-adic string solutions. None of the solutions we find, should be said, 
appears to represent tachyon matter. 

In order to explain our results we briefiy review the p-adic string model. The model 
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is defined by an action S for a scalar field. Formulated with an arbitrary spacetime 
dimensionality, the model has a parameter p, which is initially taken to be a prime number. 
The action, however, makes sense for other values of p, and even in the limit p —>■ 1, where 
it reduces to the tachyon action of [|l^] as noted in In this paper we will take p to 



be an integer that is greater than or equal to two. The p-adic action is given by 

1 , in , 1 



S = I d^xC = \ I d^x 
9p 



2 p+l 



4-4^- (1-1) 



The infinite number of spacetime derivatives are manifest in the differential operator p 2 ° . 
Here, we have 

□ = -^ + V.V. (1.2) 

and one defines 

1 °° 1 1 

p-^°=exp(--lnpn) = In p)"- . (1.3) 



n=0 



The equation of motion following from the action is rather simple looking: 

p-h°(f) = (f)P , (1.4) 

For the case of time dependent but spatially homogeneous solutions - the only kind we 
will study in this paper - the equation is 

p59?0 = 0P. (1.5) 
Finally, we note that the scalar potential is given by 

The cases odd p and even p are qualitatively different (see figure ^, although both have 
potentials unbounded from below. For odd p the potentials are even in 0, while for even p 
they fail to be even. In both cases there is a local maximum at = 1, where the tachyon 
has = —2, and a local minimum at = where the field seems naively to have no 
dynamics as its mass-squared goes to infinity. 

In order to study effectively these theories and their solutions we need expressions 
for the components of the energy-momentum tensor. For the case of a higher derivative 
theory we extend the conventional Noether procedure and find an expression for T'^ as 
an infinite series involving all numbers of derivatives. These expressions are formally local 
and can also be used for open string field theory. The energy of a solution does involve 
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Figure 1: The p-adic string potentials. For odd p the potential is even and has two 
unstable maxima. For even p there is only one unstable maximum. 

infinite number of derivatives and it is conserved. Its evaluation can sometimes be done 
in closed form. It is also determined unambiguously from the above Noether method. 
The Noether calculation of T/, where from we must read the pressure, has improvement 
ambiguities and we use the metric variation of the covariantized p-adic action to calculate 
its value. 



We also discuss the nature of the initial value problem for equations such as (|T] 
containing infinitely many time derivatives. While an equation of motion containing up 
to derivatives requires N independent initial conditions, we show that in the case 
of infinite number of time derivatives the infinite number of "initial conditions" are not 
independent but are actually subject to (possibly) infinitely many constraints. It may even 
happen that the space of consistent initial conditions is parametrized by a finite number 
of parameters. If this were the case we would have a surprising example of a theory with 
infinitely many time derivatives whose solution space is not altogether different to that 
following from equations with two time derivatives. 



The first time dependent solution in p-adic string theory was found by Brekke et. al. [0 



The equation of motion ( p..5|) was studied by expressing the differential operator in the 



left hand side as a convolution of with a gaussian, as will be reviewed in section ^ 
The equation of motion then becomes a non-linear integral equation. For the case of odd 
p the authors found numerically a kink solution of the integral equation. The solution 
interpolates in time between the vacua (f) = ±1. No solution was obtained for even p. 
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In this paper we study in detail time dependent solutions of p-adic string theory. We 
use a combination of analytic and numerical methods. In particular we show how to use 
the convolution approach to constrain solutions and to rule out certain behaviors. Some 
of the numerical methods are based on considerations of |lOl and the recent work 



While we have no exact closed forms for any solution, we have analytic expressions valid 
in certain limits. As the solutions are sometimes quite surprising we try to check them 
in two ways. We verify they hold numerically both as solutions of the nonlinear integral 
equation, and after evaluating several derivatives of the solution, we confirm that they 
directly satisfy (|1.5| ) with the derivatives expanded as in ( p..3|) . This not only increases 
our confidence on the solutions, but it is of some theoretical significance in that the 
exponentiated differential operator and the convolution form do not agree when acting 
on certain kinds of smooth functions. 



After confirming the kink solution of |T6| to high accuracy and checking it does satisfy 
both versions of the equation of motion, we turn to finding families of solutions for odd 
p potentials. We find nontrivial field oscillations around the tachyon vacuum = 0. 
These are rather surprising. It is well known that the linearized equations of motion of 
the tachyon around the vacuum admit no solutions. This is the p-adic version of the 
familiar statement that there are no conventional open string excitations at the tachyon 
vacuum |^2[. The solutions we find exist because they solve the complete nonlinear 
equations. The oscillations are completely anharmonic and the basic frequency goes to 
infinity as the amplitude of the oscillation goes to zero. The solutions are of the form 
~ A(cos{ujt)y^^, with uj"^ ~ — \n{A) and provide an explicit example of a non-standard 
open string excitation capable of carrying energy at the tachyon vacuum. These solutions 
are obtained by a level analysis of the differential equation and were checked also by 
the integral form. There are two possible interpretations for these oscillations. In the 
first one, they would represent oscillating solutions with a nontrivial energy-momentum 
tensor which would be expected to radiate quite efficiently into closed string modes, thus 
converting the tachyon energy into closed strings. In the second interpretation these 
oscillations, suitably quantized would be the closed string excitations themselves. We will 
discuss both possibilities, but believe the former is more likely to be correct. 

In the studies of tachyon matter the tachyon vacuum is at infinity and the tachyon 
never reaches this point. Moreover, there is no meaning to configurations beyond the 
tachyon vacuum. On the other hand, in string field theory or in p-adic theory, the tachyon 
vacuum is at a finite point in the field configuration space and there are field values beyond 
the tachyon vacuum. In the solutions we obtained the field typically speeds by the tachyon 
vacuum without difficulty. The case of even p, in particular p = 2, where the potential is 
cubic resembles open string field theory. For even p we have found no constant amplitude 
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oscillations around the tachyon vacuum (but we lack a proof that they cannot exist). For 
p = 2 a rolling solution down the unstable maximum at = 1 is seen to zoom by the 
tachyon vacuum = and after an excursion to negative values having potential energies 
higher than that of the original unstable vacuum, ever growing oscillations occur. In 
this solution the pressure does not go to zero for large times. The analog for even p of 
the family of odd p bounded oscillations is a family of solutions that has ever growing 
oscillations (ever growing oscillations also seem possible for odd p). 

Last but not least we study the rolling problem in the full open string field theory 
using the analytic continuation of the marginal deformation solution of as advocated 
in . We find that the tachyon rolls down towards the vacuum, goes beyond it by a large 
factor, and that a pattern of growing oscillations appears to set it. This surprising result, 
in many ways, motivated our full analysis of the p-adic string model. Indeed, the results 
summarized above concerning the p-adic model indicate that rolling in OSFT appear to 
be remarkably similar to rolling on the p = 2 potential. This lends credence to the idea 
that the OSFT solution can be trusted and is not an artifact of some approximation 
scheme. On the other hand it also suggests that this solution is not tachyon matter. 

Some concluding remarks and observations, as well as a detailed list of open questions 
can be found in the conclusion section. 



2 Convolution and Initial Value Problem 

In this section we begin by considering the convolution form of the p-adic string equation 
of motion. In this form, the equation is a nonlinear integral equation with a gaussian 
kernel. The convolution form and the differential form of the equation of motion appear 
to be equivalent only in the space of real analytic functions. 

We then show that the convolution form of the equation of motion is a useful tool 
that allows qualitative analysis of solutions. We use it to show that there cannot exist 
solutions of p-adic string theory where the field approaches monotonically the tachyon 
vacuum for large times. We also use it to find constraints on possible oscillatory solutions 
and to rule out certain types of lump solutions. 

Finally, we discuss the crucial issue of the initial value problem in a theory with infinite 
number of derivatives. 



2.1 Convolution form of the p-adic equation of motion 

In this subsection we discuss the convolution form of the p-adic string dynamics obtained 
by [|l^. Let us first derive this form. To this end we recall the differential form of the 
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equation of motion (pT5| 



p5^^?(/, = 0P. (2.1) 

As a first step we assume a well defined Fourier transform (j){k) of the p-adic field 

1 1"°° 



^{k) = ^ I e''' <j>{t)dt . 



giving (f){t) by inverse Fourier transformation: 

1 /"OO 



= I e-'''' (j){k)dk . 
2n J -co 



Then, the left hand side of the equation of motion can be written as: 

yln J—oo 



dk ^ / e~^e-''=(*-* ^dt' (j){k) 



OO 



2n J-oo \ -\/27r \np J- 

' e-—v^'-'^(l){t')dt'. (2.2) 



V^vrlnp 

Defining the gaussian convolution C[0] 

2 n T) ^ ' 



1 

■\/27r Inp 7- 



CX) 



C[0](t) = ^=^/ e-^p^'"''>^(P{t')dt' , (2.3) 



OO 



the equation of motion can then be rewritten as 

<p{tY = m{t)- (2.4) 

This is a nonlinear integral equation. Note that the convolution operator is correctly 
normalized to represent the differential operator: acting on a constant it gives the constant 
back 

C[a]=a. (2.5) 
Constant solutions (p{t) = 0o of the p-adic field equation ( |2.4|) therefore arise when : 

00 = (Po 



p>l. (2.6) 



This gives (po = and 0o = 1 valid for all p > 2. For odd p we also have 0o = —1- These 
solutions all correspond to the tachyon sitting on a critical point of the p-adic string 
potentials. 
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Before we go on and solve either or ( |2.4| ), we should note that there are smooth 
functions iplt) [i.e. continuous functions with continuous derivatives of all orders) for 
which the convolution does not represent the action of the differential operator: 

p-HV(t)^C[^](t). 

A simple example is provided by the following smooth function with compact support: 
where a and b are positive numbers such that b > a. When a < \t\ < b, i/) interpolates 




-b -a a b 



Figure 2: A smooth function that is identically one for \t\ < a, and identically zero when 
\t\ > b. For such a function the differential form of the equation of motion does not 
coincide with the convolution form. 

smoothly between and 1. Since all the derivatives of ip are zero at t = 0, we have 
p~^°^p{0) = ^/'(O) = 1; but on the other hand C['?/'](0) < 1 since C[l] = 1 and ^jj is strictly 
smaller than one when |t| > a. Therefore equation (|2.2|) does not hold for this smooth 
function. We expect, however, that ( p.2|) holds for real analytic functions whenever the 
convolution exists. We have checked many of the solutions obtained in this paper and 
verified that they solve both forms ( |2.1j ) and ( |^.4D of the equations of motion. This 
seems to mean that solutions of the p-adic string theory lie on the space of real analytic 
functions. 

2.2 Convolution constraints on solutions 

General features of a solution to the equation of motion can be gleaned from the properties 
of the gaussian convolution. We will illustrate this by proving three results that give us 
insight into the questions we are trying to address. 
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The first result rules out the possibility that the tachyon may roll monotonically down 
from 0=1 reaching the tachyon vacuum = 0. 

Claim 1: There is no solution (f){t) such that (f){t = — oo) = 1 and (f){t = +oo) = with 
(f){t) decreasing monotonically in time. 

Proof: Consider a monotonically decreasing function (j){t) satisfying the above conditions, 
as illustrated in figure ^ Since the function approaches zero, there is a time to such that 
(pito) = a, with a > and sufficiently small such that 

< I . (2.8) 

On the other hand the equation of motion requires 

a^ = Cmto). (2.9) 

Consider now the auxiliary function iplt) such that ?/'(t) = a for t < to, and ^/'(t) = for 
t > to- It is clear from the monotonicity property that (j){t) > tplt) for all t. Moreover, ip 
has a very simple gaussian convolution. Using (|2.9| ) we have 

aP = Cmo)>Cmo) = ^. (2.10) 

This violates the inequality in ( |2.8| ) and proves the claim that no solution exists satisfying 
the stated conditions. 

The second result concerns the possibility of solutions bounded in time. These could 
be bounded oscillations or otherwise. The ranges of variation are constrained by the 
following claim: 

Claim 2: Consider a solution (f){t) and constants a < b such that a < (f){t) < h for all 
times t. Moreover, assume that the values a, and h are actually attained for some specific 
times. Then < 6 < 1 and a < 0. For odd p we must also have a > —1. 

Proof: Let t^ denote a time when 0(to) = b and ti denote a time when (f){ti) = a (see 
Figure ^). From the equation of motion we must have 

6^ = C[0](to), a^ = C[0](ti). (2.11) 

Since 0(ti) < (p(t) < 4>(to) we have 

F = CM (to) < C[0(to)] =b, aP = C[0](ti) > C[0(ti)] = a, (2.12) 

giving us the inequalities 

F<b, aP>a, b>a, (2.13) 
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Figure 3: A monotonically decreasing field configuration going from the top of the po- 
tential all the way to the tachyon vacuum. Such field configuration cannot satisfy the 
equation of motion. 



where the last inequality is that from the definition of parameters. We rewrite them as 

6(y'-^-l)<0, a(a^'-^ - 1) > , h>a. (2.14) 

Note that p — 1 is a positive integer. Consider the a inequality. If both factors are positive 
we must have a > 1. If both factors are negative there are two cases: (i) even p requires 
a < 0, and (ii) odd p requires — 1 < a < 0. Consider now the h inequality. If 6 > then 
0<6<1. If6<0 there is no solution for even p, and for odd p we see that actually 
fe < — 1. The inequality h > a allows now selection of the final ranges. Since < 1 in all 
cases, the possibility a > 1 cannot be realized. Thus a < always. If p is even then a < 
and < 6 < 1. If p is odd then — 1 < a < 0. This is not consistent with h < —1, and 
therefore we get again < 6 < 1. This is the statement in the claim. 

Our third result will concern the absence of certain lump solutions for the case of even 
p. In order to establish this result we first note that for any bounded function the 
derivative of its convolution equals the convolution of its derivative: 



dC[<f) ] 



W = c5 (*). 



dt 



(2.15) 



This result follows by noting that ^ equals to ("J/) acting on the gaussian appearing 
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Figure 4: A bounded field configuration attaining a maximum value of b and a minimum 
value of a. The field equation gives constraints on the possible values of a and b. 



in the convolution (|2.3|) . This derivative can be integrated by parts without picking any 
boundary contributions and one can let it act on the function (p thus obtaining the result. 

Let us now define a monotonic lump associated to a bounded solution as in claim 2. 
This is a solution having 0(— oo) = 0(oo) = b and 0(to) = for some unique to- In 
addition decreases monotonically for t < and increases monotonically for t > t^. An 
example of such (f) is shown in figure |^. 

Claim 3: There are no monotonic lump solutions for even p. 

Proof: Because of claim 2, we have 0(to) = < 0. Since has a minimum at to, the 
convolution C[0](t) = 0^ must have a maximum at to- This requires 

-^{to)<0. (2.16) 
We differentiate (|2.15|) once more to obtain 

(to) = ^== ^ / e 2inp(* *o) (2.17) 



(it^ -\/27r Inp Inp 



oo 



dt' 



From the monotonicity condition we see that the factor in brackets is positive both for 
t' < to and for t' > to. Therefore the above integral is positive, in contradiction with 
( |2.16 ). This establishes claim 3. 
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Figure 5: A monotonic lump with endpoint values b and unique minimum value of a. This 
kind of lump solution cannot exist for even p. The dashed line represents the gaussian 
convolution of the field configuration. 

2.3 The fate of the initial value problem 

In this subsection we will make some basic observations about the initial value problem 
for the equation 

p59t0 = 0P. (2.18) 

In doing so, we view this as a differential equation for time evolution. This equation has 
an infinite number of derivatives and it is therefore unclear how to define the initial value 
problem. 

We will argue that an equation as the one above is qualitatively very different from a 
differential equation having a large but finite number of time derivatives. In solving the 
initial value problem for a differential equation having time derivatives up to order N, 
one must specify initial conditions - the values of the field and those of its first (A^ — 1) 
derivatives at an initial time, for example. Typically, those values are unconstrained, and 
solutions regular for some finite time exist for any choice of initial conditions. 

In an equation such as that in (|2.18|) one would be led to believe that one can specify 



an infinite number of independent initial conditions. More precisely, an infinite number 
minus one - that is, if we specify all derivatives of at zero time, we could use the 
equation to read the value of (j) at zero time. But a puzzle arises - if the function and all 
of its derivatives are specified at zero time, the natural assumption of analyticity would 
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imply that the function is known and can be reconstructed from its Taylor expansion. It 
seems the differential equation is not needed anymore! 

We will see below that equation ( p.l8|) actually imposes an intricate set of conditions 
for the initial conditions! We believe the analyticity assumption, and thus think that the 
evolution is determined fully by the values of the function and all of its time derivatives 
at the initial time - nevertheless finding consistent initial conditions is now the problem. 
In the system with infinite number of time derivatives the problem is not that of evolving 
an initial value configuration, but rather that of constructing a consistent initial value 
configuration. 

In an equation with up to derivatives, and with initial conditions, taking further 
derivatives of the equation typically does not yield constraints for the chosen initial values. 
The resulting equations simply determine the higher order derivatives at the initial time. 
On the other hand for equation ( p.l8| ) taking arbitrary numbers of time derivatives yields 
new constraints that must be satisfied by the infinite set of initial conditions. This is 
what makes the problem different. 

To show that there are indeed many constraints implied in ( |2.18D consider the following 



example. Let the first derivative of be non-vanishing at zero time, but all higher 
derivatives vanish: 

0'(O) = a^0, (f)"{0) = (f)"'{0) = . . . = . (2.19) 

Consider now equation ( p.l8| ) for p = 2. Consistent with the above initial conditions we 
find 

p59?0(O) = (/,(0) = 02(0) . (2.20) 

Consider now choosing 0(0) = 1 to satisfy this constraint. Since we now have the field 
and all of its derivatives specified at zero time, the (analytic) field is thus determined 
to be = 1 + at, and this does not solve equation ( |2.18[ )! What went wrong? We did 
not have a consistent set of initial conditions. Indeed, taking two derivatives of equation 
( P:TH|) we have 

p'^'^^d^m = 20'(O)0'(O) + 20(O)0"(O) , (2.21) 
which at zero time gives a new constraint on the initial conditions 

= 2a^ (2.22) 

This constraint indicates that we should not have expected a solution for nonzero a. It 
also shows that derivatives of the equation of motion yields new constraints on the initial 
values. 

It does not seem altogether unreasonable to expect that equation ( p.l8|) together with 



all of its derivatives impose an infinite number of consistency conditions on the initial 
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values. Could it be that that consistent initial values are parametrized by a finite number 
of parameters? More exploration is necessary to answer this question. 



3 Energy-Momentum in Higher Derivative Theories 

Both p-adic string theory and string field theory have actions where fields appear with in- 
finitely many derivatives, both spatial and temporal. It is also the case that the action of 
p-adic string theory has no explicit dependence on the spacetime coordinates. The same 
occurs for open string field theory in many important backgrounds. It therefore follows 
that we expect to be able to define a conserved energy-momentum tensor through a gen- 
eralized Noether procedure. Since the Noether procedure does not yield automatically a 
symmetric energy momentum tensor, we do not necessarily obtain the energy-momentum 
tensor that acts as a source of the gravitational field. This causes no comphcation for 
the case of the energy, which is unambiguous. But there are total derivative ambiguities 
in identifying the pressure from the Noether construction of {—T-). Therefore pressure 
is calculated by doing a metric variation in the relativistic covariantization of the p-adic 
action. We focus on the values of energy and pressure in the case of time dependent but 
space independent backgrounds, and then specialize further to the p-adic string model. 

3.1 Generalized Noether construction 

The situation in case is described by an action S of the type 



where the dots denote the dependence on fields acted upon by an ever increasing number 
of derivatives. The symbol in the above action may represent a collection of fields, 
but for notational simplicity we will not include an extra index. We emphasize that this 
collection of fields can include, in addition to scalars, fields of arbitrary tensor type. This 
is indeed required for the results to be useful in string field theory. All the results below 
are trivially extended when is a collection of fields. 

For a Lagrangian without explicit coordinate dependence the transformation 



where e" is a spacetime constant, gives rise to a conserved current. Under this transfor- 
mation the Lagrangian varies into a total derivative 





5(j) = e'^dacp, 



(3.2) 



SjC^e"daC^da{e"C). 



(3.3) 
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If the lagrangian C, only contained first order derivatives the stress tensor would be given 

as 

^' = -'"'^ + ^*- 

Here we are using the notation 

0Q = 9„0, (3.5) 
or more generally, for any object A, we define 

The familiar result in (|3.4| ) holds for any field of any tensor typeQ. This is because 
under constant space translations any tensor field transforms as indicated in (|3.2|) . What 
we wish to have now is the generalization of ( p. 41) for the case when the action contains 
all possible numbers of derivatives. Our result will be an expression for the stress tensor 
involving all numbers of derivatives, and reducing to the above if C only depends on first 
derivatives. 

We begin the derivation of the general formula by writing the equation of motion 
following from the action in (p.l|) . One finds 



^ dC ^ dC \ ^ dC \ dC \ 



In here the derivatives are defined as 



with no additional symmetrizations. Thus while (pf^-^^-.-fj,,, is by definition symmetric in all 
indices, tjt-^ — is not. It is also useful to note that 

dC dC dC . . 

'^^^ = ^ + 'P^^^ + ^ <Pc,viV2 H ^ ^7 (pav^-vk ^ ■ (3-9) 

This equation encodes the lack of explicit coordinate dependence in the lagrangian. 

At this stage it is simplest to state the result for the energy-momentum tensor and 
then confirm it is conserved. We find 



^For fields other than scalars the resulting energy-momentum tensor may not be symmetric. 



16 



dC 



dC_ 
dC 



+(-!)'{(; 



dC 



dC 






) (t>av^ } 
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dC ^ 





■'av2 



+ 



dC 



dC 



+ (-!)'( 



^UlU2--llUk 



4>auiiy2---iyk } 



UlU2---Uk-l 



We must now verify that 



9 = 



(3.10) 
(3.11) 



when we use the equation of motion ( |3.7|) and equation (|3.9|) . The verification is straight- 
forward once we note that each group of terms in between braces has a simple derivative. 
When taking a on any such group the only contributions are from: 9^ acting on the C 
derivative in the first term, and, from 9^ acting on the field derivative in the last term. 
All other derivatives cancel each other due to the alternating signs in the set of terms. 
As a result we find 



d 



dC 



^UlU2tl 



+{(-1)1- 



l/ll/2fJ- 



dC 



^. + G 



dC 



■'a^v\V2 



^I/lI/2---I'fcM 



uiiy2---i^kfJ- 



dC 



^p.VlV2---fk 



4>aiMUii/2---i^k} 

I (3.12) 

Using the equation of motion (|3.7|) one finds that the terms multiplying (pa add up to 
Using ( ^l9|) one then sees that indeed d^T^ = 0. 

As written, the expression for the energy-momentum tensor can be used directly to 
obtain its form in p-adic string theory. As mentioned before, this result can be readily 
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extended to a collection of fields (ff. In this case one simply sums over the various fields 
by replacing in ( |3.10D the terms 



I^^i (3.13) 



With this replacement we can use ( p.lO|) to find the energy-momentum tensor associated 
to a solution in open string field theory. 

3.2 Energy in time dependent solutions 

The focus in this paper is on solutions that have time dependence but no spatial depen- 
dence. For these class of solutions the expression ( p.lO|) for the energy momentum tensor 



simplifies considerably and all the tensor components can only have time dependence. 
From the conservation law, and with j, denoting spatial indices, 

= dX = doT^ + ^^T^ = doTl^ ■ (3-14) 

Since the energy Tq (t) is conserved we do not have ambiguities in its construction - an 
additive constant ambiguity is fixed by the condition that when the field does not vary 
in time, the energy must coincide with the potential energy. 
For solutions with no spatial dependence we also have that 



0, (3.15) 



since all space derivatives must vanish. Therefore, the terms involving field derivatives in 
( |3.10|) vanish when calculating T° and Tj. We find 



T^ = 0, Ti = -S]C. (3.16) 

Note that the conservation law d^T-" = gives doT^ + djT- = or equivalently djT^ = 0. 
This conservation does not rule out time dependent improvement terms in T/. Therefore 
the identification of the pressure as p{t) = —Tu = —T22 = ■ ■ ■ = +C{t) is ambiguous. 
This is also clear from the fact that it is set to equal the Lagrangian, a quantity for which 
total time derivatives are ambiguous. We will give a computation of the pressure using 
the methods of general relativity in the following subsection. 

Let us now calculate the energy density Tq . Denoting, for any quantity A{t) 

^n-^^' (3-17) 
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we find from ( |3.1CI| ) that the /c*'' term in braces becomes 

^-'Hi7^).i'^-i-Sr-). ■^^ + "- + (-i)'(7il^)„M' (3-18) 

and can be rewritten as 

k 

We then recognize that the complete expression for the energy is simply 

oo k 

^ = To° = -/: + EE(-ir(7r^) (3-20) 



Shifting the sum over k by one unit we write 



CO fc-l O/' 

E{t) = -C + E E(-l)"^(o-)^0.-..- (3.21) 

k=l in=0 '^Vk 

This is the result we needed. It provides an expansion for the energy in terms of time 
derivatives of the field. If the lagrangian depends only on finite number of derivatives 
then the expansion terminates for some finite value of k. If the lagrangian has infinite 
number of time derivatives then the expression for the energy is an infinite series. That 
series may or may not be possible to sum in closed form. Nevertheless the expression for 
the energy is formally local - it gives the value of E{t) in terms of derivatives all of which 
are evaluated at t. 



We will now evaluate the above expression for the energy in the case of the p-adic 
string theory ( |1.1|) . For any solution p'^i^cf) = (jf oi the equation of motion the lagrangian 
density evaluated at the solution becomes 



C=-^\\-^c^^\ (3.22) 
2 1 + p 



In addition, the derivative Mr- only exists for even k, and reads: 



|^ = -A,0(-M^1. (3.23) 



With this result, the expression for the energy in ( ^.21| ) becomes 



1 ^2^1 1 A .1 1 



Eit) = -£ + ^ - ( Inp) ^(-1)™0^ he-m, (3.24) 
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or, after minor rearrangements, 



21-1 



'->m (P2l-m ■ 



(3.25) 



m=0 



This formula will be used to confirm the correctness of certain solutions of p-adic string 
theory by testing energy conservation. We will also use it in section ^ to compute the 
energy of oscillatory solutions. 



3.3 Calculation of the pressure 

In this section we use the unambiguous definition of the energy-momentum tensor from 
general relativity. We make the p-adic model action invariant under coordinate transfor- 
mations by including a spacetime metric gap and calculate the energy-momentum tensor 
as 



T, 



5S 



The action of the p-adic model, with derivatives expanded and covariantized reads 



(3.26) 



9p 



+ 



p+l 



1 ^/ 1 nU 



. (3.27) 



Here the box operators are the covariant ones, and since all of them act on scalars they 
can be all written as 



a. 



-A, 



-99"" d,' 



(3.2^ 



The terms that are delicate to vary involve repeated action of box operators: 

1 



J dx^J~g(f) (j) = j dx (f) df^i-^- 



-99'"''d,, 



-9 



-gg^'^'d,. 



1 



-9 



(3.29) 



where all derivatives act to the right. Since we are interested in pressure we will only 
consider the computation of the components of the energy-momentum tensor Tq,/? having 
spatial indices. Since the solutions we focus on are only time dependent, this means 
that we need not vary the metric components g'^^'^^ appearing above for these would yield 
contributions to the pressure that have explicit spatial derivatives. It therefore suffices to 
vary the factors of \^—g, which is done with the help of 



6^ 



1 

'2^ 



-9 9al3 ^5'"^ 



6- 



1 



1 1 

277= 



(3.30) 
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When we vary a particular ^—g in ( p.29| ) we must integrate by parts all the derivatives 
to the left of the variation. There are (2£ — 1) such variations in ( |3.29D , each one giving a 



different splitting of the derivatives between the two fields. All in all we find (with % not 
summed) that the pressure is given by 

11 1 1001 „i 21— \ 

Pit) = -Tl = -l (-2'^^ + ^'^^"O + ^ ^ (2 ^"^) ^ • ^^-^^^ 

i/p ^ ^ ^i/p i=i ^ ra=\ 

This pressure formula will be used to analyze rolling solutions. 



4 Rolling Down the Unstable p-adic Vacuum 

In this section we construct solutions representing rolling of the tachyon from the unstable 
vacuum. That is, we want solutions where in the infinite past the tachyon approaches the 
unstable vacuum. We will obtain a solution for the case of odd p (even potentials) and 
another very different solution for the case of even p. 



It was already pointed out in ||T6[ that ( p.4| ) has a kink and an anti-kink solution when 



•p is odd. These are the solutions we are interested in for they start on one maximum 
at t —00. As it turns out they end up on the other maximum at t ^ +00. We 
will construct such solutions to high accuracy, and verify that they not only solve the 
convolution form of the field equation, but also the differential form. 

The case of even p is much more surprising. Here one could have expected a lump 



solution, but as we saw in section p.2| these are not expected to exist. Essentially the 
difficulties arise because even though the field can go into negative values the convolu- 
tion, that must equal an even power of the field, cannot. What we find (for p = 2) is 
rolling down, an overshooting and then ever-growing oscillations. This is the result of a 
calculation using an ansatz for the solution in the form of a sum of exponentials. The 
result is so surprising that we check it in two ways. One by verifying that the solution 
thus constructed satisfies the convolution form of the field equation, and second, by eval- 
uating the energy and checking that it is conserved for the indicated motion. We also 
calculate the pressure and find that it does not go to zero. We will see in section ^ that 
this behavior seems to occur also in OSFT. 



4.1 Kink solution for odd -p 

Even though the analytical form of the kinks are not known, it is relatively easy to 
construct them numerically. For example, one can use the iterative procedure where a 
given field configuration 0(t) is used to calculate the function C The field equation 
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then implies that (C[0])^/^ gives a new (and possibly improved) approximation to (j){t). 
Thus we hope for a situation where the iteration 

0(t) (^72^ / e-^(*-*')V(t')c^t') , (4.1) 

converges to some definite answer. Note that for odd p there is no complication taking 
the root - positive numbers are taken to have positive roots, and negative numbers are 
taken to have negative roots. 

Starting with the step function 



<P{t) 



+1, t<0, 

0, t = 0, (4.2) 
-1, t>0, 



(f){t) will converge to the kink solution when we iterate (|4.1| ). Figure ^ shows the kink 
solutions for p = 3, p = 5 and p = 11. As can be seen from ( |2.4| ) (and as already noted 
1^), near = the field behaves like 0(t) ~ t^^^. Also when p oo, the kink tends 



m 



to a step function. Note that this means the field zooms by the tachyon vacuum crossing 
it with infinite velocity and spending the least possible amount of time in its vicinity. 

Our numerical kink solutions have been constructed from ( |2.4|) . As mentioned before, 
it is important to check that they also satisfy equation ( |2.1|) . From equation ( p..3| ) we have 
that 

■i° = exp(ii„pa?) = i:.(^)'"^9?'"- («) 



p 



m=0 



We are thus led to define an approximate version D^^") of this differential operator where 
we include time derivatives of order less than or equal to 2n: 

D^'^'Uito) ^tC-^Y^ (dtr <P{to) . (4.4) 
^oV 2 y ml 

In Table |l| we perform a check of the solutions obtained for p = 3. We choose two values 
of the field = 0.999, 0.99 to do our check. For these values we compare 0^ to D^'^'^^cp, 
where the derivatives are evaluated using the numerical solution for the p = 3 kink. If 
0(to) satisfies (^), then D^'^''^{to) should converge to 0^(to) when n oo. In the first 
case, only eight derivatives can be evaluated accurately; in the second case, we can trust 
only six derivatives. In both cases, this suffices to see a relatively good agreement between 
D(^")0(to) and 0^(to)- We thus conclude that the kink obtained by the convolution form 
of the equation of motion satisfies also the differential form. 

If we try to repeat the same check for smaller values of the field, we see that high 
derivatives have large absolute values, and it would then be necessary to sum a large 
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-1 



Figure 6: Kink solutions for j9 = 3 (solid line), p = 5 (dotted line), and p = 11 (dashed 
line). The horizontal axis is time, and the value of the field 0(t) interpolates between the 
unstable vacuum (p = 1 in the far past and the unstable vacuum = — 1 in the far future. 



to 


0(to) 


/^^'V(to) 








0'(to) 


-4.67086 


0.999 


0.997905 


0.997309 


0.997067 


0.996974 


0.997002 


-3.04298 


0.99 


0.978765 


0.972588 


0.971867 




0.970299 



Table 1: Checking eqn. (|2.1| ) on the numerical kink solution {p = 3) computed from (|2.4|) . 

number of them to observe convergence. Finer numerical methods would be needed to 
evaluate enough derivatives. It will become clear in section ^ that the kink is the zero- 
frequency limit of solutions describing anharmonic oscillations; and this will give us further 
evidence that the differential form of the equation is indeed satisfied for the kink. 

4.2 Rolling tachyons for p-adic strings with even p 

The case of even p is rather different from the case of odd p. When p is even the potential 
has only one maximum (0=1) and we cannot have a kink solution. We could expect that 
the kink gets replaced by a lump-like solution, but as we showed in section a typical 
lump solution does not exist. We will see that if the tachyon rolls down the unstable 
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maximum towards the tachyon vacuum it overshoots it and then ever growing oscillations 
begin. We will first argue qualitatively why such behavior occurs. Then we will calculate 
the field evolution using a suitable expansion. Finally we will test the solution in two 
ways. 



4.2.1 Qualitative discussion of rolling solution 

We now imagine the tachyon starts to roll down from (f) = 1 in the infinite past towards 



the tachyon vacuum = 0. Because of claim 1 of section we know that the field 
cannot go asymptotically to the tachyon vacuum for large time. Thus we expect the field 
to go past = into negative values, as illustrated in figure |^ where at some time to 
we have 0(to) = 0. As soon as the field turns negative it becomes delicate satisfying the 
equation of motion, since the convolution must remain positive. This would be very hard 
to achieve if the field remains negative for a long time. We can therefore expect the field 
to reach a minimum at some time ti and then go back to zero value at some time t2- 
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Figure 7: A rolling solution for even p and a consistency analysis based on the value of 
the convolution C[0] shown in dashed lines. 

In figure ^ a dashed line shows the convolution, which must equal 0^. The convolution 
at to must be zero, and this is possible since the field is substantially negative to the right 
of to- But right after to the convolution must increase. How can that happen when the 
gaussian peak is moving towards values where is turning more negative ? This can 
happen if after t2 the field 0(t) is becoming positive very quickly. It must become so large 
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and positive that even the fast-decreasing gaussian tail manages to pick up a significant 
contribution. 

Now consider the convolution at ti. It can be positive because of the large positive 
values of the field beyond t2. But this convolution must start decreasing to the right of 
ti. It must do so even though the gaussian peak is moving towards values of the field that 
are large and positive. This can only happen if shortly after t2 the field turns negative 
again and does so extremely fast. All in all we find a pattern of ever growing oscillations. 
We will now confirm this somewhat qualitative analysis with a computation. 



4.2.2 Rolling with ever-growing oscillations 

We now look for a solution where the tachyon rolls down from its maximum. We shall 
make the following ansatz 

oo 

0(t) = l-^a„e""*, (4.5) 

n=l 

where a > is a constant to be determined. Note that we are using "harmonics" of the 
basic exponential e"*. This is sensible due to the structure of the equation we are trying 
to solve. With a positive, moreover, we guarantee that as t —>■ —oo the field is at the 
maximum. This kind of expansion has also been discussed in |^ . 

It is possible to anticipate the value of a. It should correspond to -\/— M^, where 
is the value of the mass-squared at the maximum. This can also be derived from the field 
equation, where we will confirm that ai plays the role of a "marginal" parameter. Since 
= — 2 at the maximum (for any p) we must have a = a/2. 

The ansatz (|4.5|) for the solution is plugged into the left hand side of (^]l|) and for 
p = 2 we find 

oo 

p59* 0(t) = 1-J2an 2^"'"' e""* . (4.6) 

n=l 

In addition, the right hand side of ( p.l| ), letting ao = —1 takes the form: 



oo n 



n=0 m=0 

Therefore the equation gives the set of relations 



(0W)' = E(E «ma„-^)e""*. (4.7) 



n 

12 2 



-a„22°" =^a^an-m, n>l. (4.8) 

m=0 

For the case n = 1 we find the condition: 

- ai25"' = -2ai , (4.9) 
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which imphes that a = \/2, and Oi is arbitrary. The relations in ( [4 .81) then give 

n-l 

an{2''' -2) = -Y. aman-m- n>2. (4.10) 

m=l 

It is clear from this equation that we can determine the exact values of the coefficients a„ 
iteratively. For example, we readily find 

a2 = -^. (4.11) 

It may seem at first sight that we have the freedom to choose Oi, but a shift in time can 
be used to rescale it to one (or to minus one if oi is negative, in which case the tachyon 
rolls towards the unbounded side of the potential). This is possible because ~ a", as 
can be seen from (|4.1CI|) and ([4.11|), and the structure of the ansatz (|4.5|). 




Figure 8: The first oscillation of 0(t) as function of time. The field oscillates with ever- 
growing amplitude. 

The above iterative procedure can be carried to any desired order to find a solution 
that is accurate for longer and longer times. Here we show the result for (j){t) to level five 

0(t) = 1 _ e^* + — e^v^* - ^ e=^^* 
^ 14 3570 

+ '^-^ e^v^* e^^* + • • • . (4.12) 

3275389320 6105767870038200 ^ ^ 
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to = 


to = 1 


to = 2 


to = 3 


to =4 


to = 5 


0(to) 


0.0711485 


-1.92423 


3.17788 


184.953 


-424.634 


-196021 




0.00506212 


3.70267 


10.0989 


34207.7 


180314 


3.8424- 10^° 


(to) 


0.00506212 


3.70267 


10.0989 


34207.7 


180314 


3.8424- 10^° 



Table 2: Comparison between 0^(to) and C [0] (to), where is the solution to level 30. 
The agreement is remarkable. 

The signs of a„ alternate, that will make the field oscillate. In figure |^, we plot the first 
oscillation. We see that the tachyon starts with the value one, rolls slowly towards the 
minimum, then overshoots the point = — | - the point where the tachyon would turn 
around in a theory with the same potential and with a canonical kinetic term. Here the 
field goes down to a value of about = — 2.6. In doing so the field has climbed a height 
which is about 55 times that of the original height of the unstable minimum! This can be 
consistent with energy conservation because in the p-adic string model the kinetic energy 
can be negative. The field turns around at = —2.6, overshooting the maximum and 
continuing to grow until it reaches a value of about = 593, a very large value indeed. 
At this point the potential is very negative and therefore the kinetic energy must be very 
large and positive. The field then manages to climb up the potential to go to lower values, 
and continues to oscillate with ever-growing amplitude. 

4.2.3 Testing and exploring the solution 

We want to do two kinds of check on this solution. First, since it was constructed from 
the derivative form of the equation of motion, we want to make sure that it satisfies the 
convolution form of the equation of motion. This will also confirm the correctness of 
the intuition that suggested that convolution requires ever growing oscillations. Then we 
want to measure the energy of the solution. In particular, we want to confirm the rather 
unusual property that ever growing oscillations can be compatible with the constancy of 
the total energy. 

In table |, we compare 0^(to) to the convolution C [0] (to) calculated numerically. We 
have taken 0(t) to level 30. We see that the agreement holds to at least six digits accuracy. 
This is impressive given that the field values include very large numbers and therefore are 
sampling at least two cycles of oscillation. 

In the second test we use the energy formula ( p.25| ) to see if the solution satisfies 
the requirement that the energy is constant. Since the energy formula includes infinite 
number of derivatives we expect better accuracy as we include more and more terms. 
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Figure 9: Testing the constancy of the energy E{t) for the rolhng solution with ever 
growing oscillations. This figure shows gpE{t), calculated from ( p.25| ), using a level 30 
solution ([4.5|). The successively flatter (solid) curves are calculated from ( p.25|) , keeping 
respectively 40, 60, 100, and 150 time derivatives. The dashed curve is 0.014>(t) and is 
just shown reference. 



This is indeed what we observe in figure We note that including more derivatives in 
the energy formula makes it constant for longer times (when the higher derivatives of 0(t) 
become more important). It is also clear on this graph, that the energy of the rolling 
solution is the height of the maximum of the potential. 

In figure |10| we show the pressure (multiplied by Qp) as a function of time. This was 
calculated using ( |3.31| ) and ( [4. 12] ) evaluating the solution to level 14, and using up to 
80 derivatives in the pressure formula. Note that for large negative times the value is 
negative and equal to (—1/6). Shown in dashed lines, for comparison is the pressure 
including only the first derivative contribution. We see no sign that the pressure is going 
to zero for large times, therefore no evidence that this solution represents tachyon matter. 

It is of interest to compare this solution with similar rolling solutions in ordinary field 
theory. As noted in [21| the perturbative expansion of a rolling solution in a 0^ theory 
using the exponentials associated to the unstable maximum breaks down at a finite time. 
This does not happen with the p-adic string solution. The expansion coefficients a„ fall 
off extremely fast due to the factor of 2"^ in ( |4.10 ). 
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Figure 10: Calculating the pressure of the maximal rolling solution for the p = 2 potential. 
Shown in dashed lines is the pressure (multiplied by g^) as calculated including only first 
derivatives in ( |3.31| ). In continuous line is the pressure including all derivatives that are 
relevant in the interval shown. 

5 Anharmonic Oscillations Around the Vacuum 

The p-adic string model, as reviewed in the introduction, satisfies one important property 
expected from the tachyon. At the minimum, the mass of the field goes to infinity and 
there are no conventional degrees of freedom. Indeed, for 0(t) small, the linearized equa- 
tion of motion is ]9~2°0 = 0, which requires □ = +oo, that is, = oo. We therefore 
have no conventional harmonic oscillations at the bottom of the potential. 

It will be seen, however, that while there are no oscillations that solve the linearized 
equations of motion, there are oscillatory solutions of the nonlinear equations of motion. 
This is the case when p is odd, and the potential is even. Such oscillations do not appear 
to exist when p is even. 

The oscillations represent a family of solutions - in the limit as the frequency goes to 
zero they go into the kink solution discussed earlier. In the anharmonic oscillations to 
be constructed below the amplitude of oscillation is a function of the frequency. As the 
frequency goes to infinity, the amplitude will go to zero. These oscillations can be found 
numerically, either by solving the convolution equation by iteration ( [4.1D starting from a 
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periodic (f)o{t), or by a "level" truncation analysis of the differential form of the equation 
of motion, applied on a Fourier expansion of 0. We will explain this last method in detail, 
because it will be useful in order to understand the analytic form of in the large 
frequency limit. It will also allow us to calculate both the amplitude/frequency relation, 
and the energy of the oscillations in the large cu limit. 

5.1 Series construction and amplitude/frequency relation 

Let us write a Fourier series for (f){t) which we imagine oscillates with some amplitude A 
and period T around the tachyon vacuum = of an odd p potential. First, we choose 
the origin of time axis such that (f)(0) — A and 9^0(0) = 0. Also, since the equation of 
motion is invariant under (f) —(f), we can demand that —(pit) = (f){t + ^), where T = 
with uj the fundamental frequency. This implies that (f> can be written in terms of odd 
modes only 

oo 

(f){t) = Yl «2n+l COS {{2n + l)ujt) . (5.1) 
n=0 

Since cu is defined to be the fundamental frequency we assume ai ^ 0. It is easily seen 
that 

oo 

P^'Uit) = E «2n+i COS ((2n + l)u;t) . (5.2) 

n=0 

The right hand side (f)^(t) of the equation of motion can also be expressed as a sum of 
cosines with odd modes only because p is odd. 

Now wc want to truncate 0(t) to level (2A^ + 1): 0(t) = J2n=o '^2n+i cos ((2n + l)ujt), 
and for concreteness we will take p = 3. To level one, (f>{t) = aiCos{ujt). Plugging this 
into the equation of motion gives 

12 3 1 

p~^^ Qi cos{u!t) — (oi cos{u!t))^ ~ 4 ^1 cos(a;t) + - cos(3a;i) , (5.3) 

and, since we work to level one, we keep only the cos(a;t) term: 

1 3 

p-^'^\i = -al. (5.4) 

This equation has the nontrivial solutions 

al = ±-|3-i-^ (5.5) 

With UJ = 1, this is ai = ±0.877383. Note that at this level ai is the amplitude of 
oscillation, and it depends nontrivially on the frequency. To level three, (f>{t) = ai cos(a;t)+ 
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as cos(3ci;t), and we have two equations 

_i££ 3 3 2 ^2 /r 

p 2 ai = -a^ + -a^as- -aia^, (5.6) 

p 2 as = -a^ + -fliag - -ag. (5.7) 

With cij = 1, we have, apart from trivial solutions and obvious sign flips, 

ai = 0.930367, ag = -0.153804 . (5.8) 

With these values the amplitude A ~ 0.777. This is not too far from the value A ~ 0.818 
obtained by a rather accurate iterative convolution calculation. In figure |TT], we show 
the solutions to level 1, 3 and 9, compared to the very accurate solution obtained from 
convolution. At level nine the two different numerical methods give almost the same 
solution. This means that the oscillations satisfy both the convolution and the differential 
form of the equations of motion. Since the — limit of the anharmonic oscillation is 
the kink, we also expect that the kink is a solution of both equations of motion, as verified 
to some degree in section |4.1| . 

We emphasize that in this example we started by choosing a frequency u and the 
numerical solution was then found to be unique (up to obvious time shifts). In particular, 
the amplitude is given uniquely in terms of u. In addition, there are solutions for any 
value of the frequency. This situation is drastically different from the case of classical 
harmonic oscillations in conventional field theory, where the frequency is given by the 
mass of the field, but the amplitude is free to take any value. 

Although we do not have a closed form solution describing the oscillations, we can 
derive such closed form in the u —>■ oo limit. Indeed, in this limit, the right hand side of 
^21) is dominated by the first term in the sum: 



p 2-°0(t) =a^p-Vcos(cut) + 0(p-5'^') . (5.9) 
From the equation of motion, we thus have 

(j)P{t) = ai cos (cut) + O (p-t'^') . (5.10) 
In the a; — cxD, the 0{..) terms can be neglected and the above equation implies that 

0(t) -_z!^ A[cos(cut)]'/P , (5.11) 
where the amplitude A is related to ai and to uj as: 

AP = aip-^ . (5.12) 
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Figure 11: Solutions for the anharmonic oscillation with p = 3 and = 1, to levels 1 
(dashed line), 3 (dotted line), 9 (thin solid line), and the more accurate solution from 
the iteration procedure (thick line). The horizontal axis is time. The amplitude is A = 
0.817615 

Recall that since p is odd, there is no problem in taking the root in ( ^.11|) - for negative 
arguments the root is defined to be negative. The evaluation of ai is possible from the 
approximate form of 0(t). Indeed ai is the first harmonic in the expansion of [cos{ujt)Y^P. 



Using ( |5.1|) we have 

A [cos(co't)]^^^ = ai cos{ujt) + cos{3ujt) + ■ ■ ■ , (5.13) 
and therefore we can calculate 

ai = — I cos(ujt)(l)(t)dt "^-^^ — / (cosx) dx 

n Jo n Jo 

= A^\-^. (5.14) 



^r(| + ^) 

Now combining ( [5.12| ) and ( p.l4| ), we find the asymptotic expression for the amplitude: 



2rfl + f) 

\2 + 2p) 
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UJ = 0.6 


UJ = 0.8 


UJ = 1 


UJ = 1.5 


uj = 2 


Numerical amplitude 


0.962253 


0.899904 


0.817615 


0.580459 


0.358948 


Amplitude from (|5.15|) 


0.975466 


0.903262 


0.818225 


0.58046 


0.358948 



Table 3: Comparison of the amplitudes of the p = 3 anharmonic oscillations, calculated 
from the numerical solutions and from the asymptotic formula (|5.15|) . 



To illustrate the asymptotic forms we show in figure |T2] one period of oscillations for 
various values of uj and the asymptotic profile ( |5.11|) . 




Figure 12: One period T of anharmonic oscillations for p = 3. On the first row from left 
to right: uj = 0.1, uj = 0.2, and uj = 0.3. On the second row from left to right: uj = 0.5, 
UJ = 2, and the last graph shows the asymptotic profile (cos cut) The scale has been 
adjusted to show all periods as equal. 



At last, in table 0, we compare the amplitudes of some p = 3 numerical solutions 
(whose absolute accuracy is about 10~^) to the asymptotic amplitudes ( |5.15| ). We see 
that both the profile and the amplitude approach rapidly their asymptotic limits. The 
asymptotic amplitude is already precise to six digits when uj >2. 
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5.2 Energy of anharmonic oscillations 

It is possible to simplify the formula ( p.25| ) for the energy in the case of anharmonic 
oscillations. We will eventually find the following asymptotic expression for the energy 

E ^uj'AP+\ (5.16) 

where hp is a p-dependent constant that will be determined in the following calculation. 

Let us plug the expansion into ( p.25| ). Since the energy is conserved, we can 

evaluate it at t = 0; at this particular time, all the odd derivatives of vanish; and the 
m-th derivative, for m even, is 

oo 

0m(O) = E a2„+i(-l)^(2n + 1)V . (5.17) 



n=0 



Thus ( p.25| ) becomes 

m even 

OO OO 

X Y: a2n+i(-l)^(2n + 1)V ^ a2,+i(-l)'-^(2fc + i)^^-™,,^.— , 

n=0 k=0 

After obvious simplifications, and rearranging the order of summation, we get 
E = -C{t = 0) - ^ E a2n+ia2.+i E( ^ 7[ E ( 

m even 

The last sum is a finite geometric series, it can thus be summed explicitly 

n = k 



2^-2 /2n + l 



E 



\2k + l) 



The energy then becomes 

,1^2/ ^^2(2^ + 1)2, uj\2n + l)\ nU 

^i/p n=0 f=0 ^• 

-| OO 1 OO 2 ^ -| 

- E «2n+ia2.+l 7— -2E(-t1^p) 7^ ((2^ + If - (2^^ + 1) 



2£ 



2^7^ 1 _ ^2n±i V fr{^ 2 ^' i\ 
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We can now do explicitly the I sums and they give simply exponentials. And since is a 
solution of the equation of motion, we have 

where A = 0(0) is the amplitude of the oscillation. The expression for the energy thus 
simplifies to 



E 



We now want to find the asymptotic behavior of ( 5.18| ) when u; — >■ oo. This is easily done 

,2 



by noting, from (|5.15|) , that p "a ~ ^, and from (|5.13|) , that Oj^+i ~ A'^. With this 



we see that the leading term in the first sum of (|5.18|) is the first term, and it is of order 
^2^p+i ^YiQ other hand the leading term in the second sum is of order A^^^. This 
is also the order of the first term in the equation. Therefore, in the limit that we are 
considering, the leading term in the energy is 

9'. 4 ,^20Fr(| + 4) 

where we have used ( p.l4| ) and ( |5.15| ). In total the final expression is 




+ (9(^1^+1)+ higher powers of A , (5.19) 



where the first correction is of order A^"*"^, without an tu^ in front of it. The next terms are 
higher powers of A, and their relative contributions are thus exponentially vanishing when 
a; — »• cxD. It is straightforward to confirm that the above result would arise by computing 
the energy using simply = ai cos(a;t). Thus in this large frequency approximation the 
first harmonic carries the leading contribution to the energy. 

In harmonic oscillations, the energy goes like A'^u'^, where A is the amplitude and u 
is the frequency. The frequency dependence of the energy for the anharmonic oscillations 
is the same, but the amplitude dependence is different - the amplitude appears with the 
same power as it appears in the potential. Of course, in the present case the amphtude 
and frequency are not independent variables. 
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In summary, we have found physical, energy- carrying excitations around the tachyon 
vacuum. If we try to interpret them in the context of string theory, they could correspond 
to non-conventional open string excitations that may radiate into closed strings. In the 
Sen conjecture, where the tachyon vacuum is supposed to be the closed string vacuum, 
physical excitations around the closed string vacuum would naturally be interpreted as 
closed strings. Could the above anharmonic oscillations represent closed strings? Perhaps, 
but there are a few complications. First, these solutions do not have analogs for even 
p potentials where any oscillation must eventually grow without limit. In some ways 
the even p case seems more closely related to bosonic string theory. In addition, we 
have calculated what would be the analog of classical open string solutions. The closed 
string states would correspond to the quantization of the above oscillations. Since the 
closed string spectrum should be coupling constant independent, the quantization must 
somehow cancel out the factor l/g"^ present in the energy of the classical solutions. For 
other discussion of closed strings and the tachyon vacuum see ||24||. 



6 Family of Solutions for Even p 

We have already found a family of solutions for the case of odd p - the anharmonic 
oscillations around the tachyon vacuum. In this section we will present a continuous 
family of solutions for even p that are, in some sense the closest analog to the anharmonic 
oscillations of the odd p theory. The oscillations, however, will be unbounded. This can 
happen even in the limit when the energy can be vanishingly small. 

In studying the maximal rolling solution for even p in section [4.2| , we saw that the 
"frequency" a in the expansion (|4.5| ) was forced to take the value ^/2. We will now 
describe a family of solutions labelled by a parameter a such that < a < ^/2. For 
simplicity we will consider only the case p = 2. Our ansatz is therefore 

oo 

(pit) = ^ a„ cosh(nat) . (6.20) 

n=0 

A solution of this form has time-reversal symmetry (f){—t) = And at t = 0, the 

tachyon field has value 0(0) = I]^o ^n- Let us set up the level truncation scheme for this 
ansatz. For convenience is it useful to rewrite the expansion as 

-| oo 
n=— oo 

with 

a_„ = (in = ttn , for n 7^ , oo = 2ao . (6.22) 
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Again, the left hand side of the equation of motion is simple 



oo 



pH'0(t) = l a„p5"^"^e""*. (6.23) 



n=— oo 



The right hand side of the equation of motion gives: 

-| oo oo 

'^'W = 7 E ( E «man-™)e""*. (6.24) 



n=— oo m= 



The last two equations imply that the equation of motion is equivalent to the set of 
equations: 

12 2 1 °° 

° ^2 ^ a,„an-m- (6.25) 

m=— oo 

The structure of these equations ensures that it suffices to consider n > 0. Moreover it 
is clear that given that the sum in the right hand side is infinite one cannot solve these 
equations exactly. A level expansion, of course, is possible. To work at level means to 
keep all equations above with < n < A^ and to keep all cig with |g| < A^. Finally, let 
us note that given a solution with expansion coefficients a„ we get another solution by 
letting 

a„^(-l)"a„. (6.26) 

To level zero the equation of motion is ao = Og, and has for solution a tachyon sitting 
either at the maximum or at the minimum of the potential. To level one, we have two 
equations for ao and ai 

«o = «o + y , (6-27) 
25°'ai = 2aoai . (6.28) 

Apart from trivial solutions, one finds 



c2 , 

ao = 2--^ 



5 

2 1 / 2 



ai = ±2-^2 1/2-2-. (6.29) 

Some interesting features can already be seen to this level. First of all, note that we have 
two branches, indexed by the sign of ai, and as could be anticipated by ( |6.26|) . If ai is 
positive, a higher level analysis indicates that so will be all higher coefficients a„, and 
the solution diverges to positive field values without oscillating. When oi is negative, the 
signs of the will alternate, making the solution oscillate with ever-increasing amplitude. 



37 



a 


1.41 


1.4 


1.3 


1.2 


1 


0.8 


9iE 


0.16082 


0.147657 


0.0592576 


0.0205954 


0.00114305 


6.59996- 10^6 



Table 4: Energy of the solution for some values of the parameter a. When a is close to 
-\/2, we have QpE |. And as a decreases, the energy goes to zero very fast. 



It is also clear from ( |6.29| ) that a has to be smaller than \/2 in order for ai to be real. 
In this regard, the maximal rolling solution of section [i.2| , having a = \f2 can be thought 
as a limit case. It cannot be obtained, however, by taking the a \/2 limit of ( |6.20| ), 
that would only give a trivial solution. 

The values of the field at t = and at this present level one approximation are (taking 
Oi to be negative) 

Xxm m = 1 , limm = ^ - i= ^ -0.207. (6.30) 

a^V2 2 v2 

Higher level analysis shows that for any a between and -\/2, the solution oscillates 
without a bound. When a y^, we have 0(0) 1, and the energy tends to 7r-2, the 
height of the potential at its maximum. When a —>■ however, 0(0) — * and the field 
becomes flatter near t = 0, in other words one has to wait for longer times before we see 
it taking large values. This is illustrated in figure |13|. Moreover the energy goes to zero 
very fast as a becomes small. Even though in these solutions the field can spend a long 
time oscillating around the tachyon vacuum their energy is necessarily very small, and 
thus they cannot represent tachyon matter. In table ^, we show a few values calculated 
from the formula ( p.25|) applied on numerical solutions at level 30. 



7 Rolling the Tachyon in Open String Field Theory 

In this section we study the rolling of the tachyon in OSFT. There are at least two ways 
we could do this analysis. In the first one, we would use level expansion to focus on a few 
low levels, write the equations of motion (containing infinitely many time derivatives) and 
solve them as in the p-adic string case. A second method was proposed in ref. (section 
2), and it uses the analytic continuation of the marginal string field theory solution of 
ref.p3|. This is the method we will follow here to examine the rolling tachyon. 

We will consider the solution representing a marginal deformation of a Dl brane 
stretched along a circle of unit radius. This marginal deformation can be used to inter- 
polate from the Dl into a DO brane, and it involves tachyon and higher field harmonics 
of the type cos(nX), where X is a coordinate along the circle. The analytic continuation 
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Figure 13: Tachyon profiles as functions of t for a = 1.4 (dashed line), a = 1.2 (dotted 
line), and a = 0.8 (solid line). The curves become flatter near the origin when a is 
decreased. 



requires replacing cos(nX) — > cosh(nt). We use the solutions of PB[, section 3, and let the 



radius go to one. In addition, we work including fields up to level four and interactions 



up to level eight, using the action given in Appendix B of [|23|. In this way we can include 
up to the second harmonic of the tachyon. The marginal direction is represented by the 
first tachyon harmonic Ti appearing in the tachyon field expansion 

T{t) = To + Ti cosh(t) + T2 cosh(2t) + ■ ■ ■ . (7.1) 

The string field equations are solved by setting a value for Ti and solving for the other 
components of the string field. We take Ti = 0.05 and the equations of motion give 

T{t) = 0.00162997 + 0.05 cosh(t) - 0.000189714 cosh(2t) . (7.2) 

It is in fact this result that motivated much of the analysis in the present paper, and a few 
comments are thus in order. First note that use of the cosh expansion implies that the 
field satisfies T(0) = 0. The other initial condition could be viewed as giving the value 
of the tachyon at t = 0. But certainly that is not the way the equation is solved. Fixing 
Ti one solves for all other harmonics To,T2, ■ ■ ■, and all of them together fix the initial 
value of the tachyon. For Ti = 0.05 we find from the above equation that T(0) = 0.05144 
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Figure 14: The tachyon profile as a function of time as it rolls down in OSFT. The initial 
conditions are T(0) = 0.05144 and T'(0) = 0. The tachyon rolls past the tachyon vacuum 
T ~ 0.55 and the first turning point is at about T ~ 1.65. 

and T'(0) = 0. The first harmonic cosh(t) reflects the blowing up of the tachyon in the 
approximation where the potential is quadratic. Note that the second harmonic comes 
about with negative sign, thus at some time the tachyon stops growing and turns around. 
We had expected naively that this turning point would be around the tachyon vacuum, 
thus signaling in this approximation that the tachyon does not cross over. But the result 
does not support this expectation. In here the tachyon overshoots the tachyon vacuum 
by a large factor. The maximum is reached for t ~ 4.88 and gives T = 1.649, which is 
about three times the value T ~ 0.55 representing the tachyon vacuum. The tachyon as a 
function of time, is shown in figure III. We have verified that the turning point T ~ 1.65 is 



quite stable under small changes of the marginal parameter Ti, supporting the conclusion 
that the tachyon generically overshoots the tachyon vacuum before turning around. 

The analogy with the p = 2 maximal rolling solution is actually quite striking. Recall 
the discussion below eqn. ( [4.12| ), where we observed that the p-adic tachyon climbs a 
height equal to about 55 times the height of the unstable vacuum. We can estimate easily 
the corresponding quantity in the above OSFT computation. Taking the tachyon potential 
to be the level zero one V = — + ^^j^T^, an overshooting to T = 1.65 implies that 
the OSFT tachyon climbs up about 56 times the height of the unstable vacuum! While 
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such close agreement is most likely a coincidence, the fact that these quantities are close 
confirms that the p-adic model seems to capture quite well OSFT features in a simpler 
context. 

What happens next? To answer this we would need the value of T3 the third tachyon 
harmonic. This, however, would require a level (9,18) computation which is a major new 



task. There are indications, however, that the next coefficient in (jT^j) would be positive, 
and that ever growing oscillations occur, just as in the case of p-adic strings with even p. 

For the value Ti = 0.05 used above, the next three fields u,v and w, all of them at 
level 2, are: 

u{t) = 0.000663831 + 9.216 x lO-^cosh(t) , 

v{t) = -0.00162516 + 0.00003999 cosh(t), (7.3) 

w{t) = 0.000301312 - 6.325 x 10"^ cosh(t) . 

With such small values for the first harmonic, these fields do not appear to change sig- 
nificantly the physics of the tachyon. This was also the case in the original marginal 
deformation problem, where level expansion was seen to operate quite well. In this level 
(4,8) solution there are additional fields at level three and at level four. We will not write 
here their values. 

We will not attempt here a computation of the pressure in this OSFT solution, but 
such computation would be of interest. With the behavior of the solution strikingly similar 
to that of the p-adic string solution for p = 2, we see no reason to expect that the pressure 
goes to zero asymptotically, but there could be surprises. A computation would settle 
this important issue. 



8 Concluding Remarks and Open Questions 

In this paper we have studied tachyon dynamics in string field theory and in p-adic string 
theory facing up to the novelties due to the infinite number of time derivatives in the field 
equations. 

Several classes of rolling solutions of p-adic string theory have been studied. We 
have not been exhaustive in showing all of them nor in trying to show that we have a 
complete set of solutions. For example, there are non-oscillatory diverging solutions for 
p = 2 related by the sign change in ( |6.26| ) to the oscillatory solutions. In addition there 
are similar solutions for p = 3 (and possibly all odd p) . These solutions also show ever 
growing oscillations, and can be modified by the sign change in (|6.26|) and by an overall 
change of sign. There are also additional solutions (both for even and odd p) arising 
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from the analytic continuation of the famihar lump solutions of p-adic string theory. At 
present, for any value of the energy we know of a finite number of solutions. 
In general terms our analysis has revealed a few surprises: 

• The solution space of the equations of motions appears to consist of real analytic 
(smooth) functions. Solutions appear to be in one-to-one correspondence with con- 
sistent initial values for the field and all of its derivatives. The dynamical equation 
imposes a possibly infinite number of conditions on the initial values and thus re- 
stricts considerably the solution space. This remarkable feature of the equation of 
motion was not anticipated by the consideration of systems with finite number of 
time derivatives. 

• There are energy-carrying solutions of the field equations around the tachyon vac- 
uum. These bounded but anharmonic oscillations shown to exist in the case of 
p-adic strings with even potentials, represent solutions of the non-linear equations 
not anticipated by the linearized (cohomology-like) problem that admits no solu- 
tions. 

• There exist oscillatory solutions with ever-growing amplitudes and constant energy. 
In fact, such solutions exist even for energies that approach zero. These solutions 
exhibit behaviors not seen in lagrangians quadratic in first time derivatives where 
the kinetic energy is positive definite. In the systems considered here the kinetic 
energy can be negative and thus one can see the tachyon move to higher and higher 
heights on the tachyon potential while conserving the total energy. 

At the technical level we have calculated the energy-momentum tensor in higher deriva- 
tive theories, finding useful expressions for the energy and the pressure in the case of 
time-dependent solutions. We have also learned how to extract qualitative behavior of 
the solutions using the convolution form of the field equation. Finally, we have shown 
how to obtain solutions by numerical methods, and by analytic methods in certain limits. 

Interesting open questions remain. We list below a few of them: 

• From the physical viewpoint, the most puzzling result has been that none of the 
rolling solutions obtained here appear to represent tachyon matter. That is, we have 
not found solutions where with varying values of the energy the pressure goes to 
zero for long times. The rolling solutions described in this paper would seem to yield 
oscillating energy-momentum tensors that could effectively convert the D-brane en- 
ergy into closed string excitations. This had been the conventional wisdom (see, for 
example |^) before studies of tachyon matter. Given the exciting possibility that 
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a stable form of tachyon matter could have astrophysical consequences it seems of 
utmost importance to confirm its (theoretical!) existence using string field theory. 
Since the tachyon vacuum is at a finite and apparently regular value in the field 
space of string field theory, tachyon matter must be, if present, an exotic solution. 

For the case of even p-adic potentials, the anharmonic excitations around the 
tachyon vacuum do represent non-conventional type of excitations. The most straight- 
forward interpretation is that they are degrees of freedom that would radiate their 
energy into closed strings. Could they actually represent closed strings? While pre- 
liminary indications do not support this interpretation, further investigation of the 
physics and interpretation of these solutions are of interest. 

We have seen that there is no conventional initial value problem in the sense that 
finding consistent initial conditions appears to be the same problem as finding the 
solution. The space of consistent initial conditions appears to be strongly con- 
strained in the nonlinear equations of motion of p-adic string theory. How many 
parameters does the space of consistent initial conditions have? 

Our work has given an explicit construction for the energy density of spatially 
homogeneous solutions in OSFT. The calculation of the pressure was more subtle 
and was only done for the p-adic string. The OSFT pressure should be computed 
to test if, as opposed to the case of the p-adic model, the ever growing oscillations 
could have an asymptotically vanishing pressure. A more geometrical construction 
would use open/closed string field theory where the energy-momentum tensor 



would arise from a variation of the closed string field. It may be of interest to give 
such stringy construction. 

One feature of string field theory is that the infinitely many time derivatives of 
the covariant formulation turn into first order time derivatives in the light cone 
formulation. Studies of gauge fixing in string field theory have not revealed how this 
transformation actually happens. It would therefore be of interest to understand 
if there is a light-cone gauge p-adic string theory. If the answer is affirmative the 
passage from covariant to light-cone formulation could be tractable and provide 
much insight. 

The convolution form of the field equation seems to hide completely the causality 
properties of the time evolution. It would be of interest to formulate and explore a 
test of causality in p-adic string theory. 
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It is clear that in p-adic string theory and in OSFT we have a rich set of problems and 
issues that should allow us to improve our understanding of time evolution in string theory. 
Such understanding seems essential for the development of string theory cosmology. 
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